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We examine the conformations of a model for short closed 
DNA. The molecule is represented as a cylindrically symmet- 
ric elastic string with a constraint corresponding to a specifi- 
cation of the linking number. We obtain analytic expressions 
for configuration, elastic energy, twist and linking number for 
a family of solutions spanning from a circle to a 'figure eight'. 
We suggest ways to use our construction to make other con- 
figurations and models relevant to studies of DNA structure. 
We estimate the effect of thermodynamic fluctuations. 



I. INTRODUCTION 

The elastic model of DNA has been the subject of 
intense research in the past 30 years. The approaches 
taken include Lagrangian mechanics Jl|-^| , b-splines and 
(numerical) molecular dynamics ||, and statistical me- 
chanics 0. Until now no one has been able to find the 
general equilibrium solutions for a bent twisted segment 
of the molecule. As a particular case the conformations 
of closed circular DNA have remained elusive. 

In this paper we present a method for obtaining the 
shape, twist, elastic energy and linking number of an 
isotropic elastic segment subject to constraints. We will 
use this formalism to describe unknoted circular DNA. 
We will also outline how our formalism can be exploited 
to produce other shapes of interest to DNA researchers. 



II. ELASTIC MODEL AND EXPRESSIONS OF 
INTEREST 

In a previous work we outlined some developments in 
the elastic model of DNA Q]. The molecule is repre- 
sented as a slender cylindrical elastic rod. At each point 
s we describe the rod by relating the local coordinate 
frame C to the frame Cq rigidly embedded in the curve 
in its relaxed configuration. The relationship between 
the stressed and unstressed local frames is specified by 
Euler angles 0(s), cj)(s),ip(s) needed to rotate Cq into C. 
Let us summarize some of the results needed for this pa- 
per. Let the elastic constants of bending and torsional 
stiffness be denoted, respectively, by A and C, and let L 
be the lengh of the rod. The elastic energy is given by 
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A major relevant result from our previous work is that 
Ck, the linking number, can be written as 
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Where we have used Fuller's theorem M to obtain the 
writhe, and White's theorem for the total link. We 
shall use these local quantities to determine the be- 
haviour of the family of closed solutions that extremize 
the elastic energy. 



III. CANDIDATES: THE ELASTIC EXTREMA 

A. A Family of Curves 

We are searching for a family of solutions which are 
extrema of elastic energy of a closed molecule. This fam- 
ily of 'writhing' solutions will be indexed by the writhe, 
Wr, of a member - ranging from Wr = for a sim- 
ple circle to Wr — 1 for a 'figure eight'. We know that 
for Ck < \/3A/C the molecule remains stable as a cir- 
cle ]i|JL2[]. Therefore one bounding member of the fam- 
ily must be a simple circle, say, the XY plane with a 
Tw — Ck = ^jZAjC. We shall see that on the other end 
the bounding member, a 'figure eight', must have Tw = 
in order to be an extremum of elastic energy. However, 
the reader must remember that the actual constraint im- 
posed on the closed curve is that of a fixed Ck. We shall 
therefore examine the conditions under which the family 
of writhing curves can satisfy an imposition of a linking 
number constraint. 



B. Euler-Lagrange Minimization 

We want to extremize elastic energy and keep the curve 
closed. We find that the functional to consider for closed 
configurations of DNA is 



H 



dE P 
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(4) 
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The first term of (fl]) is clear - the elastic energy must 
be extremized. The second term will enforce the closure 
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in the z direction only because we can always find a ref- 
erence frame where the curve is closed in the XY plane. 
Written explicitly as a function of Euler angles (0) be- 
comes 



H = 



A 



ds ^ [<f sin 2 e + e 2 
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(5) 



The extrema are found by applying Euler-Lagrange equa- 
tions to (|B|). Denoting the conserved quantities as p^ = 
|J and p^sgwe get 



P4, - Pip cos 9 
A sin 9 
cos 9 



(6) 
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The equation for 9 is a quadrature obtained by integrat- 
ing §jj = ^ ^jj with Eq as the constant of integration. 
Defining u = cos 9 we see that the behavior of solutions 
is governed by a cubic polynomial in u: 
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where we order the roots c < u(= cos 9) < b < a 

(|J) requires cos 9 to oscillate between c and b. For a 
circular solution the tangent t will go back and forth 
twice. 




u = cos 9 

FIG. 1. The tangent oscillates c— >fe— >c^b^c. The 
curve is composed of two symmetric parts. 

All the relevant quantities, including the shape of the 
curve, can be obtained using (0). For example: 
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A. Constraints 

To complete the analysis we must determine the pa- 
rameters generated by the minimization procedure: the 
invariants p^ andp^,, the constant Eq, and the Lagrange 
multiplier F. To determine the parameters we impose 
constraints on the curve. Figure [j] is a good visual guide 
to the geometric meaning of the constraints that follow: 

J* ds = L/2 (fixing the length at L) (10) 

J 2 u(s)ds = (enforcing closure in z) (11) 

L 

Jo </>{s)ds = 7r ((f) must go all the way around) (12) 



B. Reparametrization 

At this point we reparametrize the problem by F, a, 6, c 
instead of F,Eq,p^,,p^. Parametrizing the problem by 
the roots of the polynomial is extremely advantageous: 
it makes the analytic manipulations more transparent; it 
also streamlines the computational tasks. The two sets 
of parameters are related in the following manner: 



E 



P<t> = 



F(a + b + c)) 



(13) 
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— (pi±p 2 ); Pi> — Y \Pi T P2) (14) 

with p/is = [(c ± 1)(6 ± l)(a ± 1)] 1/2 

The choice of branch (±) is imposed by the family of 
configurations sought. For circular DNA without intrin- 
sic curvature, takes a — , and p^, correspondingly a +.F] 
Let us make some definitions that make notation more 
transparent: 
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Employing (p|) we rewrite the constraint equations, (|l 
( |TT| ) , and ( |l2| ) , in terms of the new parameters 

32A 
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(17) 
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Where K, E and II are complete elliptic integrals of the 
first, second and third kind, respectively. 



IV. CONSTRAINTS AND PARAMETERS 

the choice of a particular branch is a non-trivial procedure 
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C. Solving for a, b, c and F 



(y/a + 1 - Van - l) = 1 



(21) 



We find that the optimal procedure for calculation of 
solutions' parameters is to use ( |l7j ) to eliminate F, then 
( |l8| ) and ( fl9|) eliminate Aa and c. Because the Wr is a 
monotonically increasing function of Ab, we leave it to 
classify the family of curves. Once the constants F, a, b, c 
are determined, the desired solutions and all the relevant 
quantities are computed via elliptic integrals. For exam- 
ple, the explicit expression for 9(s) in the first quarter of 
oscillation is (we have inverted 



cos 9 = Ab sn 2 



'FAc 
~2A l 




(20) 



(j) and ip are obtained similarly from (^|) and (f^). 
The following figure shows the family of curves computed 
in the manner discussed above. Since the constraint 
equations involve elliptic integrals, finding a solution on 
a computer is virtually instanteneous.^ 





Mr 1 

FIG. 2. The family of curves ranges from the circle in 
the XY plane to the 'figure eight' in the Y Z plane. The 
perspective is slightly assymetric to aid visualisation. 



D. Bounding Members: Circle and 'Figure Eight' 

Let us check whether the initial member of our family, 
a circle with Wr = joins smoothly with the previously 
known stable family of twisted circles (l^Jl^]. The circle 
corresponds to Wr = 0. A circle in the XY plane the 
curve must have 6q = c o = 0. ([n]) now states that 



which in turn gives 

2nA 1 , . . . 2ttA r- 

P*o = -r (v^T+T + vV^T) = (22) 

Li V a 

Combining (||) and (R) to obtain the Tw of the circle p 
gives 



Tw 
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(23) 



It is also of interest to compute the Tw of the 'figure 
eight' which, like the circle, can be done virtually by in- 
spection. The figure lies in the YZ plane, which forces <j> 
to behave as follows: (refer to Figure [I] for visualization) 



5(0) + 6 



(24) 



= which 



Combining (|24| ) and (||) forces p<p a — p^, & 
immediately sets Tw s = 0. The value of A6 8 is easily 
determined from the fact that Cs = — 1 (this can be seen 
from the curve itself) necessitating a — 1. Then jl9| ) 
gives A6 8 = 1.6522.... 



V. LINKING NUMBER AND THE 
PLECTONEMIC TRANSITION 

We have found a family of writhing solutions that are 
the extrema of elastic energy. The writhe of the curves 
covers the range < Wr < 1. [] However, the physical 
constraint imposed on the molecule is Ck, the linking 
number. The explicit expressions for Tw, Ck and Wr 
are easily obtained: 
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and, using White's theorem [10] 
Ck = Tw + Wr 



(25) 

(26) 
(27) 



Plotting Ck vs. Wr for our family we can discern three 
distinct types of behavior depending on the value of the 
ratio A/C. 



2 analytically and computationally elliptic integrals are 
equivalent to, say, arcsin 



3 That the 'figure eight' has Wr = 1 just before crossing can 
be seen just from the shape of the curve. 
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A/C » 1 




FIG. 3. The three types of behaviour of Ck for our 
writhing family of curves. The solid line indicates conforma- 
tions realized by the molecule, whereas the dashed regions are 
forbidden by the imposition of Ck. The ratio A/C controlls 
the plectonemic transition. 



of simple trigonometric functions. Q One complication 
does arise in cases of initial curvature: the solutions have 
to be minimized wrt to the initial value of ip. An intu- 
itive example of this caveat is the fact that the energy 
of intrinsically curved loops is changed by 'rolling', un- 
like that of the configurations which are straight in their 
relaxed state. 



B. Trefoils and Other Torus Knots 

Another interesting generalization of our results is 
obtained by altering the closure condition |]l2| , |l6| . 
For example, a curious configuration is attained by 
requiring the curve to close after <j> has turned 
4-7T rather than the 2tt required for the closure of a cir- 
cular conformation. The resulting curve is a trefoil knot. 
These curious shapes are probably the analytic analogs 
of similar configurations observed by Shlick and Olscn 
jfj . cite myself unpublished 



We can now describe what happens to the molecule as 
the linking number is increased. Initially the molecule 
remains a circle until ACk — y/ZA/C. After that there 
are three possibilities. If A/C <C 1 then the writhing 
family can support a steady increase to Ck = 1 of the 
'figure eight' and the molecule folds continuosly until self- 
crossing. (Futher increase of Ck presumably results in a 
plectonemic configuration). If A/C ^S> 1 then none of the 
members of the writhing family can support the necce- 
sary linking number, and as soon as ACk exceeds y/3A/C 
the twisted circle snaps into a plectoneme. An interesting 
case is the intermediate one, where, at first, the molecule 
starts writhing, but at some point, unable to support fur- 
ther increase in ACk it snaps into a plectoneme. Another 
interesting fact is that this behavior is independent of the 
length of the molecule. 



VI. POSSIBLE GENERALIZATIONS 

A. Intrinsic Curvature 

The method of construction of the circular DNA can 
be easily generalized to other configurations of inter- 
est. Considerable activity has been generated around 
the studies of DNA with intrinsic curvature. We feel 
such generalizations are important because various chem- 
ical effects - hydrophobic, etc. - can be modeled by in- 
trinsic curvature |T3|~|l5|] . Inclusion of intrinsic curvature 
into the present framework involves a simple modifica- 
tion of the LHS of equation jl9| ) [Q . If the inital curva- 
ture is such that the DNA is fully relaxed in the circular 
state, the characteristic polynomial (@) degenerates into a 
quadratic and all the quantities are expressible in terms 



VII. ENTROPIC CORRECTIONS 

Although our model is applicable to the cases when 
thermodanamic fluctuations are not important (i.e. short 
segments) it is worth while to estimate the contribution 
of thermal fluctuations to the free energy. Since the aim 
of this paper is to derive the elastic shapes of DNA, we 
will concentrate on thermal effects that influence curva- 
ture. This can be done in analogy to Marko and Siggia 
p7| who, in turn, used the scaling ideas originally devel- 
oped for the study of fluctuating membranes JlS| ]. Their 
(slightly generalized) dimensional analysis states that the 
entropic contribution to the free energy can be written 
as 



Af = k B T 



1/3 



k 2 = 



\A/k B T / 
where the curvature k is 
9 2 + 6 2 sin0 2 



(28) 



(29) 



We conceed that equation (28) gives only the form of the 
thermal correction, sans the multiplicative factor of order 
unity. ( p8| ) reflects the fact that increasing the curvature 
reduces the correlation length, thus suppressing fluctua- 
tion and incurring a corresponding loss of entropy. 



In the current formalism ©will force a — > inf; but it is 
easier to omit the term F cos 8 in the initial functional and 
derive the equations anew. 
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VIII. CONCLUSION 



We have presented a formalism for obtaining the elas- 
tic mimima of circular DNA subject to a constraint in the 
linking number. Aside from the simplicity of the state- 
ment of the problem, and the beauty (in the eyes of the 
creators) of the solution, our formalism can significantly 
improve the investigations of linear molecules with elas- 
tic models. This method can be easily generalized to cre- 
ate a truly 'finite' finite-element-analysis in which elastic 
models can be traced out by finite segments. The confor- 
mations of these segments can be specified in close form 
which would decrease computation time significantly. 
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